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Abstract 

Exactly solvable rationally-extended radial oscillator potentials, whose 
wavefunctions can be expressed in terms of Laguerre-type exceptional orthog- 
onal polynomials, are constructed in the framework of fcth-order supersym- 
metric quantum mechanics, with special emphasis on k = 2. It is shown that 
for \i = 1, 2, and 3, there exist exactly fi distinct potentials of fith type and 
associated families of exceptional orthogonal polynomials, where /i denotes 
the degree of the polynomial arising in the denominator of the potentials. 
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1 Introduction 



In recent years, one of the most interesting developments in quantum mechanics 
has been the construction of new exactly solvable potentials connected with the ap- 
pearance of families of exceptional orthogonal polynomials (EOP) in mathematical 
physics. In contrast with families of classical orthogonal polynomials, which start 
with a constant, the exceptional ones start with some polynomial of degree m > 1, 
but still form complete sets with respect to some positive-definite measure. 

The first examples of EOP, the so-called Laguerre- and Jacobi-type X\ families, 
corresponding to m = 1, were proposed in the context of Sturm-Liouville theory 
[U [2] . They were then applied to quantum mechanics and proved to be related to 
some exactly solvable rational extensions of known potentials, with the additional 
interesting property of being translationally shape invariant [3j H] (although at that 
time the list of additive shape-invariant potentials was thought to be complete [5J 

El EI). 

Shortly thereafter, the first examples of X 2 EOP and of related shape-invariant 
potentials were proposed in the framework of conventional supersymmetric quantum 
mechanics (SUSYQM), suggesting a growing complexity with m |8J. Two distinct 
families of Laguerre- and Jacobi-type X m families were then constructed for arbi- 
trary large m [HI ED] and their properties were studied [TT], Q2] . Several approaches 
were considered in connection with the Darboux-Crum transformation [T3l [TH [T5] , 
the Darboux-Backlund one [HJ HZ] or the prepotential method [TSJ [TH]. Possible 
applications to position-dependent mass systems [20], to PT-symmetric potentials 
[U [21], to conditionally-exactly solvable ones [221 [23], and to the Dirac and Fokker- 
Planck equations [21] were studied. 

Very recently, the X m EOP were generalized to multi-indexed families by means 
of multi-step Darboux algebraic transformations [25] or the Crum-Adler mechanism 
[26]. 

The purpose of the present work is to show how the related exactly solvable 
potentials arise in an (essentially equivalent) simple extension to higher order of the 
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SUSYQM approach used in [U |Sj. For simplicity's sake, we shall restrict ourselves 
here to rational extensions of the radial oscillator connected with Laguerre-type 
EOP. 



2 Second-order SUSYQM 

Going from first- to second-order SUSYQM may be achieved in two slightly different 
ways, corresponding to the parasupersymmetric (PSUSY) scheme [27] or to the 
second-derivative (SSUSY) setting [281 US EH EE E2] . 

In the former approach, one starts with a pair of first-order SUSYQM partners 
(in units h = 2m = 1) 

#(+) = A^A = — ^— + U (+) (x) - E, fl-H =AA* = —^ + (x) - E, 

dx 2 dx 2 ^ 1) 

A ] = — p + Mi), A = ^- + W(x), V {± \x) = W 2 (x) tW'(x) + E, 
dx dx 

which intertwine with the first-order differential operators A and A^ as AH^ + ' = 

H^A and A^H^ = H^A^. Here W(x) is the superpotential, which can be 

expressed as W(x) = —4>'(x)/(j>{x) in terms of a (nodeless) seed solution <p{x) of the 

initial Schrodinger equation 

<h +V {+ \x)] <f>{x) = E<f>{x), (2.2) 



dx 2 

E is the factorization energy, assumed smaller than or equal to the ground-state 
energy E { +) of V^, and a prime denotes a derivative with respect to x. We shall only 
consider here the case where E < E$ + \ in which occurrence 4>(x) is nonnormalizable, 
and we shall assume that the same holds true for <p~ l (x). Then and turn 
out to be isospectral [5]. 

Next, we consider a second pair of first-order SUSYQM partners and 
with the same characteristics and distinguish all related quantities from those per- 
taining to the first pair by tildes. If we choose V <y+ \x) such that V^ + \x) = V^~\x), 
then both pairs of SUSYQM partners (#(+), #H) and (H i+ \ ffH) can be glued 
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together, so that we get a second-order PSUSY system or, equivalently, a reducible 
SSUSY one. 

The latter is described in terms of two Hamiltonians 

= ~ + W(x), = ~ + W(x) : (2.3) 

dx z dx z 

which intertwine with some second-order differential operators 

d 2 d d 2 d 

A^ = — -2p(x)— + q(x), A= — + 2p(x)— + 2p'(x) + q(x), (2.4) 
dx z dx dx z dx 

as Ah^ = h^A and A^h^ = h^A\ so that the functions p(x), q(x) and the 
potentials V^ l,2 \x) are constrained by the relations 



q(x) = —p' + p 



II / I \ * 2 
2 P I P \ ° 



2p ' \2pJ 16p 2 

II / /\2 2 ^ > 

v^ 2 \x)=t2 P '+p 2 +^-(^ + c 



2p \2p J I6p 



where c is some integration constant. 



£ 
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The relation between both approaches follows from the equations = + 
h (2) = jj(-) _ ^ A \ = ^t^t ; anc i ^4 = AA. It turns out that h {1) and /i (2) are both 
partners of some intermediate Hamiltonian h = H^~' + | = — |, that the 

constant c is related to the two factorization energies through c = E — E, and 
that the function p(x) can be expressed in terms of the two superpotentials as 
p( x ) = \(W + W). 

Instead of and 0, we may start from two seed solutions 0i and 02 of the 
initial Schrodinger equation f)2.2p with respective energies E\ and E2 (less than 
Eq +S) ) and such that 0J -1 and (p^ 1 are nonnormalizable. Then, on choosing = 0i 
and = A(f>2 = W(0i, 02)/0i, so that E = E\ and = -E2, we get 

, v ^(01,02) (gl ~ ^2)0102 , . 

PW " 2W(0 1? 2 ) 2W(0 1; 2 ) ' 1 ] 

where W(0i,02) denotes the Wronskian of 4>i(x) and 02(x). Since, from (12. 5p . the 
SSUSY partner potential can be written as V^ 2 \x) = V^(x) + 4p'(a;), it is clear 
that it can be completely determined from the knowledge of this Wronskian. 



3 Radial Harmonic Oscillator in Second-order 
SUSYQM 

Let us consider a radial oscillator potential 



V l (x)= 1 -co^ + l -V±^, (3.1) 

4r Ob 



where u and / denote the oscillator frequency and the angular momentum quantum 
number, respectively, and the range of x is the half-line < x < oo. As well 
known, the corresponding Schrodinger equation has an infinite number of bound- 
state wavef unctions, which, up to some normalization factor, can be written as 

^ oc x l+1 e-^ x2 L^\lojx 2 ) oc m (z)L^(z), v = 0, 1, 2, . . . , (3.2) 

with 

z = \ujx 2 , a = 1 + 1, rn(z) = z^ 2a+1) e^\ (3.3) 

and Ly\z) some Laguerre polynomial. The associated bound-state energies are 
given by 

E$p = uj(2v + / + §) = u(2v + a + 1). (3.4) 

Motivated by the experience gained in the first-order SUSYQM approach (see 
Eqs. (2.12), (2.16), and (2.18) of Ref. 7), as well as by subsequent developments 
HOI H3] , we may consider two different types of seed solutions <p(x) with properties 
as required in Sec. 2, namely 

</>LO) = x]{z)Lt\-z) oc ^eXlf^-kc'), (3.5) 



f(z)L^(z) oc x-'e-^hi! h \\ujx 2 ), (3.6) 



with 



xl(z) = z iv°+%i; xf(z) = z-X*-%-*; (3.7) 

and corresponding energies 

E} m = -co(a + 2m + 1), Ef m = -u(a - 2m - 1), (3.8) 



respectively. Such seed solutions are related to the two families LI and L2 of 
Laguerre-type X m EOP. Note that for type II, a must be greater than m. 

Our purpose is to construct some rationally-extended radial oscillator potentials 
^,ext( a; ) with a given I by using two seed solutions 0i an d 02, as explained in Sec. 2. 
Taking into account that the order of 0i and 2 is irrelevant, there are three types of 
possibilities for the pair (0 1? <f> 2 ) and it turns out that in each case, we have to start 
from a potential V^ + \x) = Vi>(x) with some different I'. The corresponding Wron- 
skian W(0i(x), 02 (x)) can be written in terms of some //th-degree polynomial in z, 
gn(z), itself expressible in terms of a Wronskian W(f(z), g(z)) of some appropriate 
functions of z, as follows: 

(i) V i+) = V^ 2 , 01 = 0j-2,m 1 > 02 = 0L 2 ,m 2) < 7Jl l < m 2 , 

W(0i,0 2 ) =wx{-x!i_ 2 ) 2 g ll (z), 

g,(z) = W{Lt; 2 \-z),Lt; 2 \-z)), / i = m 1 + m 2 -l; (3.9) 

(ii) \/ (+) = Vi +2 , 0i = 0j+ 2 , mi , 02 = 0" 2 , m2 , 0<m!<m2<a + 2, 
W(0i,0 2 ) = u;x(x"2)V(^)> 

^(^) = yy(LLT- 2) W,^T" 2) W), ^ = m l + m 2 -\ ] (3.10) 

(iii) \/W = Vl, 1 = 0j mi , 02 = 0^, 0<mi, 0<m 2 <a, 
M0i,0 2 ) = ^x"^), 

= z w(LSSH,Lyw) - (* + «)££>(-*)4; a) (*), 

/i = mi + m 2 + 1. (3-H) 

In all three cases, we can write (provided g^ does not have any zero on the 
half-line) 

V (1) =^-i( J E 1 + E 2 ), 

= V l -uj\?k + ±z 
\ 9» 

where a dot denotes a derivative with respect to z and C = —2u, 2u, or in case 
(i), (ii), or (in), respectively. In the associated PSUSY approach, which we cannot 
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g» 



g» 
g fl 



-(E 1 + E 2 ) + C, 



(3.12) 



detail here due to space restrictions, the intermediate potential is some VJ_i jex t, 
V5+i,ext, or VJ + i iCxt potential. Note that it reduces to some bare radial oscillator 
potential in the special case where mi = 0. Furthermore, if we adopt the reverse 
order for the 0's, e.g. 4>i — ^p m2 ' ^ 2 = ^n»i lYl case (^)> the ^ na ^ potential will be 
the same, but the intermediate one will be different, e.g. some V5_i iCxt potential in 
case (in). 

Both V {1) and V {2) have the same bound-state energy spectrum, given by 

{lo(2v + 21 + mi + m 2 — 1) in case (i), 
w(2f + 21 - mi - m 2 + 5) in case (u), (3.13) 
u;(2z/ + 2Z + mi — m 2 + 2) in case (m), 

where i/ = 0, 1, 2, .... The bound-state wavefunctions ipv(x) of V"( 2 ) can be 
obtained either by acting with A on those of ip^\x) oc r)i>(z)Lv*\z), or by 
directly inserting the expression 

ti 3 \x)ocjQy n (z), n = fi + u, u = 0,1,2,..., (3.14) 

in the Schrodinger equation for V^ 2 '(x). As a result, we obtain the following differ- 
ential equation for y n (z), 

z -f-2 + ( a + 1 - z - 2z ~) T- + ( z - «)- + z ~] Vn( z ) = Q*-n)Vn(z). (3-15) 
dz 2 V S'm/ dz 9n 57* J 

(2) 

The orthonormality and completeness of ip u (x), v = 0, 1, 2, . . . , on the half-line 
imply that for any n — \i + v, v — 0, 1, 2, . . . , the differential equation (I3.15P 
admits a rath-degree polynomial solution and that when v runs over 0, 1, 2, . . . , 
such polynomials form an orthogonal and complete set with respect to the positive- 
definite measure z a e~ z g~ 2 dz. We shall denote these polynomials by L^ mi n (z), 

^Sm 1) m 2 ,n(«)» and ^Xn.m^nW in cases (*)> (")> and ( iil ')> respectively. 

It is worth noting that in cases (z) and (zz), the differential equation ( 13.151) can 
be rewritten in a slightly different form. The definitions of g^z) in (13.91) and (13.101) . 
combined with Laguerre equation, indeed lead to 

Z g = j 2 ^ ~ ( a + z )9v + V9n in case («), ^ 
1 2zg M + (a + z)^, - in case (z*z), 



where = W(Lmi 2 \-z),L^ 2 2) (-z)) or = W{L m * 2 \z),L m ? 2) (z)), respec- 
tively. As a result, we get 



d 2 ( n 9 A d g„ 
z—r + I a + 1 — z — 2z— 2a— + 2z— 

dz 2 V 9 it J dz g^ g^ 



TP- (z) 

a, mi ,rri2,n V / 



n -^a,mi,m 2 ,n( Z ) (3-17) 



and 



^T^ + ( a + 1 - 2 - 2z— - — I (g^ + g^) 

dz 2 V. 9nJ dz g^ 



1-11,11 / X 

a, mi ,rri2,n \ ) 



= {2^-n)L l l^ m2 ^z). (3.18) 
The latter equation coincides with that obtained for Ln' mi ' (z) in Ref. 24. 

4 Some Simple Examples 

In cases (i) and (ii), the lowest-degree example for g^ corresponds to mi = 0, 
m 2 = 1, leading to go = 1 or g = —1 and giving back Laguerre polynomials: 
La l fl,i, n — ^afl,i, n = ^ Q On assuming mi = 0, m 2 = m + 1 (m > 1), we get 
either # m = Lm _1) (- Z ) or 9m = -L^ a ~ l) (z), so that L^ 0m+ln = L^ m)Tl and 
^a'om+in = ^amni hence the associated extended potentials coincide with those 
already obtained from first-order SUSYQM. The next values mi = 1, m 2 = 2 do 
not provide any new result either because g 2 = L 2 a l \z) or g 2 = —l^ 2 l \—z) and 
therefore L l ^ 12 n = L l ^ 2 n and L^f 2 n = L J a2r ^ Going to mi = 1, m 2 — 3 at last 
gives rise to a new result (the same in the two cases), because g^ = [z 3 + 3az 2 + 
3(a - l)(a + 1)^ + (a - l)a(a + l)]/3 differs from both 4 a-1) (-z) and L { ~ a ~ 1] (z) 
leading to £„ 3n and L^ 3n , respectively. 

In case (in), considering either mi > 0, m 2 = or mi = 0, m 2 > produces an 
already known result L^ miA „ = ^L, mi +i,„ or L^J m2]n = £" m2+ i jn , since g mi+1 = 
— (mi + l)L^~^l(—z) or g m2+ i = (m 2 + l)L^°^i\z) as a consequence of known 



a These equalities and the following ones are of course dependent on the normalization chosen 
for the EOP. We assume here that the latter can be appropriately adjusted. 

b It is worth observing here that the last reduction has not been noted in Ref. 24. 
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relations among Laguerre polynomials. The first example not derivable from first- 
order SUSYQM is related to mi = iri2 = 1, but the resulting g% coincides with that 
obtained for cases (i) and (ii). 

We conclude that SSUSY does not lead to any new extended potential of linear 
nor quadratic type, but gives rise to a new cubic one. Its explicit form can be 
obtained from for jj, = 3, which we rewrite as 



\ /.<••:. V, ■■ I',,,,. V Ur&t = -lo{2^ + 4z 

93 



93 I 93 " " 



93 \93 



(4.1) 



by dropping the additional constant. On combining Eq. (13. 3p with the expression 
of gs(z) given above, we get 

V, (x)-^ + ^- (42) 
ViM*)- D{x) + D 2 {x y 

where 

N^x) = 12lo[lo 2 x 4 - (21 + l) 2 + 28], 

N 2 (x) = -288w[3(2/ + l)wV + 4(2/ - 1)(2Z + 3)ux 2 

(4.3) 

+ (2/-l)(2/ + l)(2/ + 3)], 

D(x) = (tux 2 + 21 + i f - 4(3cox 2 + 21 + 1). 
This result may be compared with the two cubic-type extended potentials com- 
ing from first-order SUSYQM, which are given by Eq. (14.11) with either gz(z) = 
L^~ l \—z) or g-z(z) = L% a X \z), and which can be written in the form (14.21) with 
N^x) = \2oj\oj 2 x 4 - (21 - 9) (21 + 5)], 

N 2 (x) = -144w(2/ + 5)[(2/ + 9)tu 2 x 4 + 2(21 + 3)(2l + 5)tux 2 

+ (21 + l)(2/ + 3)(2/ + 5)], 
D(x) = (cox 2 + 21 + 5) 3 - 2(2/ + 5)(3cox 2 + 61 + 11), 



(4.4) 



or 



iVx(x) = 12u[u 2 x 4 - (21 - 3) (21 + 11)], 

N 2 (x) = U4u(2l - 3)[(2Z - 9)tu 2 x 4 + 2(21 - 3)(2Z - l)tux 2 

+ (21 -3) (2/-1) (2/ + 1)], 
D(x) = (cox 2 + 21- 3) 3 + 2(2/ - 3)(3ujx 2 + 61- 5), 

9 



(4.5) 



respectively. Note that I is restricted to / > in (14. 3 p and (14. 4p . and to I > 2 in (j4.5|) . 
All three extended potentials have the same spectrum (13. 4 p than the conventional 
radial harmonic oscillator potential (13. ip . 

5 Final Comments 

What has been done in detail in second-order SUSYQM can, in principle, be gen- 
eralized to higher order k [2B1 ESJ GEH ED The construction of the EOP and 
of the related extended potentials will then be governed by the choice of k seed 
solutions 4>i(x), 02 (#) , 4>k{ x ) of the initial Schrodinger equation and by their 
corresponding Wronskian W(0i, 02, • • • , 0fc)- 

For k — 3, for instance, it is easy to see that the lowest-degree corresponds to 
H = 3, which may arise in the pure cases obtained from = Vis, 0; = 0j_ 3jmj 
or = VJ +3 , 0j = 0[+ 3jmi , with z = 1, 2, 3, mi < m 2 < m 3 , and pi = mi + m 2 + 
m 3 -3|] Since g 3 = W(L^' 3) (-z), L^~ 3) (-z), L^' 3) (-z)) = -4" a_1) («) and 
(73 = W(L^ 3) (z),L^ Q " 3) (^,4~ a - 3) (^)) = it turns out, however, 

that we shall obtain the extended potentials associated with L^ 3n (z) and L^^^z), 
respectively. 

We conclude that considering /cth-order SUSYQM with k = 1, 2, 3, . . . leads 
to exactly \i distinct extended radial oscillator potentials and corresponding EOP 
families of /ith type for /i = 1, 2, and 3. Whether this result may be generalized to 
higher values of \i is an interesting open conjecture. 

Other important points for future study are the construction of extended po- 
tentials of Morse or Coulomb type, also connected with Laguerre-type EOP (see, 
e.g., Refs. 16 and 18), as well as that of extended potentials related to Jacobi-type 
EOP. A direct proof of the shape invariance of the new potentials remains to be 
given. The existence of different intermediate Hamiltonians, as observed in the 
PSUSY approach of Sec. 3, is also worth analyzing along the lines of type A A^-fold 
supersymmetry J3TJ [33] . 

c Here we assume m\ > 0, because mi = would lead to g M 's already found from SSUSY. 
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